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ABSTRACT 



We present stellar and gaseous kinematics of the inner « 350 pc radius of the Seyfert 
galaxy Mrk 1066 derived from J and K/ bands data obtained with the Gemini's Near-Infrared 
Integral Field Spectrograph (NIFS) at a spatial resolution of ^35 pc. The stellar velocity field 
is dominated by rotation in the galaxy plane but shows an S- shape distortion along the galaxy 
minor axis which seems to be due to an oval structure seen in an optical continuum im- 
age. Along this oval, between 170 and 280 pc from the nucleus we find a partial ring of 
low <x* (^50kms _1 ) attributed to an intermediate age stellar population. The velocity dis- 
persion of the stellar bulge (<x* « 90kms _1 ) implies a super-massive black hole mass of 
« 5.4xlO 6 M . From measurements of the emission-line fluxes and profiles ([Pn]/U.1886yum, 
[Fe n]/U .2570 /mi, Pa/3 and H2/I2. 121 8 yum), we have constructed maps for the gas centroid ve- 
locity, velocity dispersion, as well as channel maps. The velocity fields for all emission lines 
are dominated by a similar rotation pattern to that observed for the stars, but are distorted by 
the presence of two structures: (i) a compact rotating disc with radius r « 70 pc; (ii) outflows 
along the radio jet which is oriented approximately along the galaxy major axis. The compact 
rotating disc is more conspicuous in the H2 emitting gas, which presents the smallest cr values 
(< 70kms _1 ) and most clear rotation pattern, supporting a location in the galaxy plane. We 
estimate a gas mass for the disc of ~ 10 7 M©. The H2 kinematics further suggests that the 
nuclear disc is being fed by gas coming from the outer regions. The outflow is more conspic- 
uous in the [Fen] emitting gas, which presents the highest cr values (up to 150 km s" 1 ) and 
the highest blue and redshifts of up to 500 km s" 1 , while the highest stellar rotation velocity 
is only « 130 km s" 1 . We estimate a mass-outflow rate in ionized gas of « 6 x 10" 2 M o yr _1 . 
The derived kinematics for the emitting gas is similar to that observed in previous studies 
supporting that the H2 is a tracer of the AGN feeding and the [Fe 11] of its feedback. 

Key words: galaxies: individual (Mrk 1066) - galaxies: Seyfert - galaxies: ISM - galaxies: 
kinematics and dynamics - infrared: galaxies 



1 INTRODUCTION 

This work is a continuing study of the central region of 
Mrk 1066 using Gemini's Near-infrared Integral Field Spec- 
trograph (NIFS) observations i n the J and Ki bands. In 
Riff el, Storchi-Bergmann & Nagar (2010, hereafter Paper I), we 
presented maps for the emission-line flux distribution and ratios 
and discussed the gas excitation and extinction, as well as the origin 
of the nuclear continuum. The main results of Paper I are: (i) The 
line emission in the near-IR lines is extended over the whole NIFS 
field (« 700 pc) and most elongated at position angle PA= 135/3 15°, 
showing a good correlation with the optical [O m] line and radio 
continuum emission. [Cavm] /I2.3220yi/m and [S ix] /ll. 2524 //m 
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coronal lines are the exceptions, being unresolved at the nucleus; 
(ii) The [Fen] XX. 2570 //m/Pa£ vs H 2 ^2.1218yL/m/Bry diagnostic 
diagram is dominated by values typical of active galaxies; (iii) 
The reddening map obtained from the Pa/?/Bry line ratio presents 
a S-shape structure with E(B - V) reaching a value of 1.7 along 
PA^ 135/3 15°; (iv) From line ratio maps we conclude that the main 
excitation mechanism of the H2 and [Fe 11] emission lines is heat- 
ing by X-rays from the central AGN. Correlations between the ra- 
dio and line-emission maps - stronger for [Fe 11] than for the H 2 - 
suggest that shocks due to the radio jet play a role in the gas excita- 
tion; (v) The nucleus contains an unresolved infrared source whose 
continuum is well reproduced by emission fro m dust with temper - 
ature -810 K and mass ~ 1.4 x 1O- 2 M . In iRiffel et all bOlCh 
we presented stellar population synthesis, which shows that within 
~160pc from the nucleus an old stellar population dominates the 
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Figure 1. Sample fits of the stellar kinematics of the nuclear region of 
Mrk 1066 using pPXF. Top left-hand panel: fit of the nuclear spectrum; top 
right-hand panel: fit of a spectrum at 075 north-west from the nucleus; bot- 
tom left-hand panel: fit of a spectrum at 075 south-east from the nucleus; 
bottom right-hand panel: fit of a spectrum at 075 north-east from the nu- 
cleus. The observed spectra is shown in black and the resulting fit in red. 



2 OBSERVATIONS AND DATA REDUCTION 

The J and Ki bands Integral Field Unit (IFU) spectroscopic data of 
Mrk 1066 were obtained with Gemini NIFS (iMcGregor et al.l2003h 
operating with the Gemini North Adaptive Optics system ALTAIR 
in September 2008 under the programme GN-2008B-Q-30. NIFS 
has a square field of view of « 370 x 370, divided into 29 slices 
with an angular sampling of 07 103x07 042, optimized for use with 
ALTAIR. 

The J band observations were centred at 1.25 yt/m, covering 
a spectral region from 1.15 /mi to 1.36//m at a spectral resolu- 
tion of « 1.7 A (from the full width at half maximum - FWHM 
- of arc lamp lines) and at spatial resolution of 0713+0702, as ob- 
tained from the FWHM of the spatial profile of the telluric star. The 
Ki observations, centred at 2.3 iim covered the spectral range from 
2.10//m to 2.26 yt/m at spectral and spatial resolutions of « 3.3 A 
and 0715+0703, respectively. The total exposure time at each band 
was 4800 s. 

The final J and Ki data cubes, obtained by combining the indi- 
vidual exposures (which included some dithering), contain ^4200 
spectra at an angular sampling of 0705 x 0705, covering the central 
700x700 pc 2 of Mrk 1066. For a detailed description of the obser- 
vations and data reduction procedures see Paper I. 



near-IR continuum (age ~10 10 yr), while in a partial ring surround- 
ing this region, the continuum is dominated by an intermediate age 
stellar population (10 8 < age < 10 9 yr). 

Mrk 1066 is an SB0 galaxy located at a distance of 48.6 Mpc 
which harbours a Seyfert 2 nucleufl for which 1" corresponds to 
235 pc at the galaxy. Previous optical long-slit spectroscopy shows 
that the high-excitation gas (traced by the [Om] emission) and 
low-excitation gas (traced by [Nn], [Sn], [Oi] and H recombina- 
tion lines) present distinct kinematics, with the former being more 
disturbed with respect to the rotation curve observed in the low- 
excitation emission lines teower etaDll995h . The authors suggest 
that the high-excitation gas kinematics is due to an outflow driven 
by the radio jet, which is co- spatial wit h the extended line emis- 
sion, along positio n angle PA= 135/3 15° dUlvestad & Wilson! 19891 : 
Nagar et al. 1991). A distinct kinematics for the high- and low- ex- 
citation gas is also supported b y optical integral field spectroscopy 
of the inner 1 kpc of Mrk 1066 JStoklasova et alll2009h . In the near- 
IR, the emission l ines present asym metric profiles revealed by long- 
slit observations dKnop et al.ll200ll) . These asymmetries have been 
interpreted has being originated in the same outflowing gas com- 
ponent observed in [Om], while most of the near-IR line emission 
(traced by the peak of the line profi les) may be due to emission 
from gas located in the galaxy plane ( Kno p et aTll200lh . 

In the present paper, we use the observations described in Pa- 
per I to obtain two-dimensional maps for the gaseous and stellar 
kinematics of the inner « 350 pc radius of Mrk 1066, which al- 
lowed us to put additional constraints on the physical scenario for 
the circumnuclear region of this galaxy. This paper is organized as 
follows: In Section 2 we present the description of the observations 
and data reduction, Sec. 3 presents the results for the gaseous and 
stellar kinematics, which are discussed in Sec. 4. Sec. 5 presents 
the conclusions of this work. 



3 RESULTS 

As shown in Paper I, in the J band we observe the following 
emission lines: [P n] XV. 14713 and 1.18861 //m, [Fen] A 1.25702, 
1.27912, 1.29462, 1.29812, 1.32092 and 1.32814/mi, Hi Pa/?, the 
He ii line at 1.16296 iim and the [S ix] coronal line at 1.25235 iim. 

In the Ki band we observe the H2 lines at 2.12183, 
2.15420, 2.22344, 2.24776, 2.40847, 2.41367, 2.42180, 2.43697 
and 2.45485 yum, the Hi Bry, the He iXl. 14999 fjm and the 
[Cavm] A2. 32204 /im coronal line. Besides the emission lines, the 
CO stellar absorption band heads around 2.3 11m are also present in 
the spectra an have been used to derive the stellar kinematics. 



3.1 Stellar kinematics 

In order to obtain the stellar line-of- sight velocity distributions 
(LOSVD) we fitted the 12 CO and 13 CO stellar absorption band 
heads around 2.3 //m in the AT-ba nd spectra using the p e nalize d 
Pixel-Fitting (pPXF) method of ICappellari & Emselleml J2004h . 
This method finds the best fit to a galaxy spectrum by convolving 
stellar template spectra wit h a given LOSVD [L(v)], repre s ented by 
Gauss-Hermite series (e.g. Ivan der Marel & Franxlll993l : iGerhardl 
Il993h : 



L(v) = 



-(1/2)// 



r„ V2^ 



^ h m Mm(y) 



(1) 



1 NASA/IPAC Extragalactic 
http ://ned w w w. ipac . caltech . edu > 
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where y = (v - V*)/cr, V* is the stellar centroid velocity, <x* 
is the velocity dispersion, v = c\nA and c is the speed of 
light. The H m are the Hermite polynomials, hm * are the Gauss- 
Hermite moments (ICappellari & Emsellem 20041). The pPXF rou- 
tine outputs the V*, cr* and the high er-order Gauss-He rmite 
moments and h^. As discussed in iRiffel et al ] d2008h and 
Winge, Riffel & Storchi-Bergmannl d2009h . the use of a library of 
stellar templates, instead of a single stellar spectrum, is fundamen- 
tal for a reliable derivation of the stellar kinematics. In this work we 
have used template spectra from the Gemini library of late spec- 
tral type stars observed with the Gemini Near-Infrared Spectro- 
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Figure 2. Stellar kinematic maps obtained from the pPXF fit. Top: centroid velocity (left) and velocity dispersion (right) maps. Bottom: h^* (left) and 
(right) Gauss -Hermite moments. The mean uncertainties are 6 kms" 1 for centroid velocity, 8 kms -1 for <x*, and 0.03 for h^* and h^. The continuous black 
line shows the orientation of the line of nodes obtained from the modeling of the velocity field. The thick contours overlaid to the velocity field, velocity 
dispersion and maps are from the HST optical continuum image obtained with the WFPC2 through the filter F606W. 



graph (GNIRS) IFU and NIFS JWinge. Riffel & Storchi-Bergmannl 
H). 

In Figure[T] we show some fits (in red) compared to the ob- 
served spectra (in black). This figure illustrate the fits at four posi- 
tions: the nucleus (top-left panel), 075 north-west (top-right panel), 
075 south-east (bottom-left panel) and 075 north-east (bottom-right 
panel). High signal-to-noise (S/N) ratio (> 20) is required to ob- 
tain a reliable fit of the observed spectra. As the S/N decreases with 
distance from the nucleus, we have replaced the flux at each spatial 
pixel by the average flux of the 9 nearest pixels, which increased the 
S/N ratio and allowed us to measure the stellar kinematics at most 
locations, except for a few pixels close to the borders of the IFU 
field. This procedure can be justified by the fact that our data cube 
is somewhat oversampled (the pixels correspond to 07 '05 x Of! 05), 
and this average is thus equivalent to use an effective aperture of 
the order of the PSF FWHM (« (7/15). 



In Figure[2 we present two-dimensional maps of the stellar 
kinematics. The central cross marks the position of the nucleus, 
defined as the locus of the peak of the Ki band continuum. The 
top left panel of Fig.[2] shows the stellar velocity field, from which 
we subtracted the systemic velocity, obtained from the modeling 
of this velocity field (see Sec. 14.11 ). The corresponding heliocentric 
value is V s = 3587 + 9 kms -1 . The continuous black line shows 
the orientation of the line of nodes - *F = 128° + 6° measured 
from north to east, and obtained also from the modeling mentioned 
above. The mean uncertainty in the centroid velocity is « 6 km s -1 . 
The velocity field shows blueshifts to the south-east and redshifts 
to the north-west and seems to be dominated by rotation, with an 
amplitude of ^120 kms -1 and kinematical centre coincident with 
the position of the continuum peak, within the uncertainties. 

In the top-right panel of Fig. [2 we present the stellar velocity 
dispersion map, which has a mean uncertainty of 8 kms -1 and 
presents values ranging from 30 to 95 kms -1 , with the highest 
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ones observed at the nucleus. A partial oval ring of low cr* val- 
ues (« 50 kms -1 ) is observed surrounding the nucleus at « 1" 
from it. The bottom panels show the higher order Gauss-Hermite 
moments h 3 * (left) and h 4 * (right). These moments measure devi- 
ations of the line profile from a Gaussian: the parameter h 3 * mea- 
sures asymm etric deviations and the h 4 * m easures symmetric de- 
viations (e.g. Ivan der Marel & Franxfl993h . The values /z 3 * and h 4 * 
vary from -0.2 to 0.2 with mean uncertainty of 0.03. The highest 
values of /z 3 * are observed from the nucleus to « 173 south-east and 
at « \" north-west from the nucleus, while the lowest values are ob- 
served over most of the north and north-east side of the galaxy. The 
fi4* map presents the lowest values at the nucleus, while the high- 
est values are observed in a ring with radius « 1" approximately 
coincident with the ring of low cr* values. 



3.2 Gaseous kinematics 

The most commonly used method to measure the gaseous kine- 
matics is the fit of the emission-line profiles by Gaussian curves. 
Nevertheless, the emission lines in Mrk 1066 present asymmetric 
profiles at many locations, which cannot be well represented by 
Gaussian curves. We have thus used G auss-Hermite series (e.g. 
Ivan der Marel & Franx| [l993 ; Riffel 2010), using the same formal- 
ism used to derive the stellar kinematics (see Section [3J} in order 
to obtain a better fit to the emission-line profiles and a more reli- 
able measurement of the centroid velocity (from the centroid wave- 
length of the profile), velocity dispersion (cr), as well as to obtain 
the higher order h 3 and h 4 mome nts. Each em ission line profile was 
fitted by the following equation teiffefcOld) : 

GH = [1 + h 3 H 3 (w) + h 4 H 4 (w)] , (2) 

cr 

where w = and a(w) - -j=e~ w2/2 , A is the amplitude of 
the Gauss-Hermite series, A c is the centroid wavelength, h 3 and h 4 
are the Gauss-Hermite moments and H 3 (w) and H 4 (w) are Hermite 
polynomials. 

We added a linear component to Eq.|2]in order to fit the con- 
tinuum under each emission-line and thus the resulting equation 
contains seven free parameters (A, A c , cr, h 3 h 4 plus two parame- 
ters for the linear equation), which have been determined by fitting 
the line profiles. The fitting of the emission-line profile was done 
using the mpfit0 routine, written in IDL0 programing language us- 
ing the Levenberg-Marquardt least- squares method, in which ini- 
tial guesses are given for the free parameters. More details about 
the fitting procedure and its implementation can be found in Riffel 

boioh . 

We have chosen the [Pn] A1.1886/im, [Fen] Al. 2570 /mi, Pa/? 
and H2 A2. 1218//m emission lines to map the kinematics of the cen- 
tral region of Mrk 1066. These lines trace the kinematics of distinct 
gaseous species (ionized gas forbidden lines, ionized gas permitted 
lines and molecular gas). These particular lines have been chosen 
in order to minimize the measurement uncertainties because they 
have the highest S/N ratios among their species. In Fig|3]we show 
a sample of fits obtained for each line profile at 075 north-west from 
the nucleus within an aperture of 0725x0725. The solid line shows 
the observed profile, the dashed line shows the fit and the doted 



2 The MPFIT routine can be obtained from the Markwardt IDL Library at 
http://cow.physics.wisc.edu/~craigm/idl/idl.html 
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Figure 3. Sample fits (dashed lines) of the observed emission-line pro- 
files (solid lines) by Gauss-Hermite series at 075 north-west from the nu- 
cleus for an aperture of 0725x0725. The dotted lines show the residu- 
als of the fit plus an arbitrary constant (for visualization purposes only). 
Top:[Pn] /U.1886//m (left) and [Fen] /U.2570//m (right) emission lines. 
Bottom: Pa/3 and H2 i2.1218//m emission lines. 

line shows the residuals (plus an arbitrary constant for visualiza- 
tion purposes). As shown in Fig. [3] the line profiles are well fitted 
by the Gauss-Hermite series and thus allow reliable determinations 
of the centroid velocity, velocity dispersion, h 3 and h 4 . 

Figure [4] shows the velocity fields obtained from the 
Gauss-Hermite fit of the profiles for the [P 11] Al. 1886 11m, 
[Fen] /il.2570yi/m, Pa/? and H 2 /i2.1218//m emission lines, af- 
ter subtraction of the systemic velocity of the galaxy (V s = 
3587 kms -1 ). The crosses mark the position of the nucleus and 
the dashed lines show the orientation of the line of nodes derived 
from the modeling of the stellar velocity field (see Sec. 14.11 ). The 
gas velocity fields are similar to each other and are dominated by 
rotation, with the south-east side approaching and the north-west 
side receding from us. Although similar to the stellar velocity field, 
the gas velocity fields are more disturbed. The most conspicuous 
feature is more clearly observed in the H 2 velocity field: a gra- 
dient from ^125 kms -1 at «075 north-west from the nucleus to 
« -55 km s -1 at «075 south-east from it along the PA^150°, which 
seems to be a compact rotating disc immersed in the larger disc 
(stellar and gaseous). 

In Figure [5] we present the gas velocity dispersion (cr) maps. 
In order to compare our cr maps with the radio structure we have 
overlaid th e contours (blue lin es) of the 3.6 cm radio continuum 
image from lNagar ctaJ dl999h on the [Fe 11] cr map. The [Fe 11] and 
[Pn] cr maps are similar and show the highest values of up to 
^150 km s -1 in regions at and around the radio structure, while the 
lowest values (« 50 kms -1 ) are observed predominantly to south- 
east of the nucleus. The Pa/? cr map presents smaller values than 
those of [Fen] and [Pn], with the highest values of ^100 kms -1 
observed approximately along the minor axis of the galaxy. The 
lowest values (<60kms -1 ) are observed in two structures resem- 
bling spiral arms, more clearly observed in the H 2 cr map. This map 
also shows that these arms - one coming from the east and north- 
east and the other coming from the west and south-west - seem to 
curve into the major axis towards the centre into a structure which 
we have identified with a compact nuclear disc in Sec. 14.21 
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Figure 4. Centroid velocity fields for the [P n] Al. 1886 /mi (top-left panel), [Fe n] i 1.2570 //m (top right), Pa/3 (bottom left) and H 2 A2.1218/im (bottom right) 
emitting gas. The mean uncertainties are smaller than 7 kms" 1 for all emission lines. The thick black contours overlaid to the [Fen] velocity field are from a 
radio continuum image from lNagar etafl Jl999h . The cross marks the position of the Ki band continuum peak and the dashed line shows the orientation of the 
line of nodes obtained from the modeling of the stellar velocity field. 



In Figures [6] and |7] we present the h 3 and h 4 Gauss-Hermite 
moment maps. For [Pn]and [Fen], the h 3 maps show negative val- 
ues to north-west of the nucleus, while the Pa/? presents negative h 3 
values at most locations (except to east and south of the nucleus). 
However the most negative values reaching « -0.3 are observed 
predominately at ~ 1" north-west of the nucleus at the edge of the 
radio jet, indicating the presence of a blue wing in the emission-line 
profiles at this location. Some regions with positive /* 3 , with values 
of ^0.1 and reaching 0.2 at a few locations, are observed to the 
north, north-east, east and south-east of the nucleus for [Fen]and 
H2, indicating the presence of a red wing in the emission-line pro- 
files at these locations. 

The J14 maps present values close to zero at and around the 
nucleus for all emission lines - indicating no significant deviations 
from a Gaussian in the emission-line profiles. At ~075 from the 
nucleus approximately along the major axis, the /i 4 maps for the 
ionized gas show regions with values ^0.15, which indicate that 



the profiles are narrower and have broader wings than Gaussians at 
these locations. 

Coronal lines are also present in the spectra of Mrk 1066, such 
as [Six] X\. 25235 /mi and [Cavm] A2.32204//m. The [Cavm] and 
[S ix] ions have ionization potentials of 12 7.2 eV and 328.8 eV, re- 
spectively JStorchi-Bergmann et al.ll2009b and their emission lines 
can be used to map the kinematics of the high ionization emitting 
gas. Nevertheless, as discussed in Paper I, their flux distributions 
are unresolved by our observations. The [Cavm] line is affected by 
the 12 CO(3-l)/i2.323//m absorption and we have thus subtracted 
the stellar population (obtained from the fit of the stellar kinemat- 
ics) from the galaxy spectra before fitting the [Cavm] emission- 
line profile. The fit gave centroid velocity V « -12 kms -1 , which 
is consistent with the systemic velocity of the galaxy within the 
uncertainties, and velocity dispersion of <r ^145 kms -1 . The fit 
of the [Six] /1 1.25702 yt/m emission line gives V ^290 kms -1 and 
<t ^70 km s -1 , but the uncertainties are large because the [S ix] line 
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Figure 5. Velocity dispersion maps for the [Pn] /U.1886//m (top-left panel), [Fen] /U.2570//m (top right), Pa/3 (bottom left) and H2 i2.1218//m (bottom 
right) emission lines. The mean uncertainties are smaller than 10 kms -1 for all emission lines. The cross marks the position of the K4 band continuum peak, 
the arrows identify the regions of asymmetry in the flux distribution, and blue contours overlaid on the [Fe 11] cr map are from the radio continuum map of 
iNaear et alJ<1999h . 



profile is contaminated by both the [Fe 11] X\ .2570 yt/m emission line 
and underlying stellar absorptions. 

3.3 Channel maps 

In order to better sample the gas kinematics over the whole ve- 
locity distribution, including the highest velocities in the wings of 
the line profiles, we "sliced" each profile into a sequence of ve- 
locity channels. In Figure [8] we present the channel maps for the 
[Fen] /ll. 2570 //m emission line for a velocity bin of ^75 kms -1 , 
corresponding to three spectral pixels. Each panel presents the flux 
distribution in logarithmic units integrated within the velocity bin 
and centred at the velocity shown in the top-left corner of the panel 
(relative to the systemic velocity of the galaxy). The central cross 
marks the position of the nucleus and the green conto urs overlaid 
to som e panels are from the 3.6 cm radio image of iNagar et al J 
(1999). The highest blueshifts of up to -500 km s are observed to 



north-west of the nucleus at ^0785 from it, approximately coinci- 
dent with a hot spot observed in the radio map. When the blueshifts 
reach smaller absolute values, of « -300 kms -1 , the flux distri- 
bution presents emission to both sides of the nucleus along the 
PA^ 135/3 15°, which is approximately the PA of both the radio jet 
and major axis of the galaxy. The highest redshifts, of up to of 
^450 km s -1 , have origin in gas located at both sides of the nucleus 
at distances of ^0785 to north-west and « 1" to south-east. 

Figure [9] presents the velocity channels for the Pa/? emission 
line. The highest Pa/? velocities - blueshifts of -350 kms -1 and 
redshifts of 320 km s -1 - are smaller than those observed for [Fen]. 
At the highest blueshifts the emission is dominated by the nucleus, 
showing fainter emission from the north-west and south-east sides, 
while most of the redshifts are due to emission from the north-west 
side of the galaxy. From blueshifts of « -100 kms -1 to redshifts 
of « 100 km s -1 the peak of the emission shifts from the south-east 
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Figure 6. /Z3 Gauss-Hermite moment maps for [Pn], [Fen], Pa/5 and H2 emission lines. The white contours overlaid to the Pa/5 cr map are from the radio 
continuum image of Na gar et all dl999h . 
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Figure 7. Same as Fig.[6lfor the Gauss-Hermite moment. 
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Figure 8. Channel maps along the [Fen] Al. 2570 /mi emission-line profile. Each panel shows the intensity distribution at the velocity shown in the top of the 
panel, integrated within a velocity bin of ^75 kms -1 . The cross marks the position of the nucleus and the green contours overlaid to some panels are from the 
3.6 cm radio image. 



to the north-west of the nucleus along the galaxy major axis, sug- 
gesting rotation. 

The velocity channel maps for the H2 /I2.1218//m emission 
line are shown in Figure [TUl for velocity bins of ~60 kms -1 , cor- 
responding to two spectral pixels. The pattern here is somewhat 
distinct from those observed for [Fe 11] and Pa/?. The H2 emission 
is more uniformly distributed over the whole IFU field, being less 
"collimated" along the radio axis, with blueshifts observed mostly 
to the south-east and redshifts to the north-west. This is what is ex- 
pected if the velocity field is dominated by rotation similar to that 
of the stars. Nevertheless, it can be noticed that the velocity fields 
are not symmetric relative to the galaxy major axis (what should be 
the case for pure rotation). For example, in the blueshifted chan- 



nels centred at -104 kms -1 and -42 kms -1 the flux distribution 
curves towards the north, while in the redshifted channels centred 
at 81 kms -1 and 142 kms -1 the flux distribution suggest some cur- 
vature to the south. These regions are marked with arrows in Fig- 
ureEOl 



The [Pn] /U.1886//m velocity-channel maps are similar to 
those for [Fen], but are much noisier and are thus not shown here. 
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Figure 9. Same as Fig.[8]for the Pa/3 emission-line profile. 



4 DISCUSSION 



4.1 Stellar kinematics 



As can be observed in Fig. [2] the stellar velocity field presents a 
rotation pattern, and in order to obtain the systemic velocity, ori- 
entation of the line of nodes and an estimate for the bulge mass of 
Mrk 1066 we have fitted the stellar velocities with a model of cir- 
cular orbits in a plane subject to a Plummer gravitational potential, 
given by 



O = 



GM 



(3) 



+ a z 



where a is a scale length, r is the distance from the nucleus in the 
plane of the galaxy, M is the mass inside r and G is the Newton's 
gravitational constant. Although this model is very simple, it has 
approximately reproduced the stellar velocities in the central re- 
gions of other Seyfert galaxies in previous studies dBarbosa et al.l 
l200dlRiffeletail2008h. 

Following Barbosa et al ] j2006h . defining the coordinates of 
the kinematical centre as (X , Y ) relative to the photometric nu- 
cleus, the observed centroid velocity at position (R, X F), where R is 
the projected distance from the nucleus in the plane of the sky and 
*F is the corresponding position angle, can be expressed as 



V r = V s + 



R 2 GM 



sin(/)cos(¥ - ¥o) 



^ 2+A2 ) 3/2 (cos 2 (^- 



¥o)- 



sin 2 or-y ) y 

cos 2 (/) / 



3/4 



(4) 



where V s is the systemic velocity, i is the inclination of the disc 
(i = for a face on disc) and is the position angle of the line 
of nodes. The relations between r and R, and between a and A are: 

r - aR and a = aA, where a = ^Jcos 2 ^ - *Fo) + " "co^^ • 

The equation above contains six free parameters, including the 
kinematical centre, which can be determined by fitting the model 
to the observed velocity field. This was done using a Levenberg- 
Marquardt least- squares fitting algorithm, in which initial guesses 
are given for the free parameters. In the fitting we considered a 
typical uncertainty in the observed velocity of 6 km s _1 . 

The parameters derived from the fit are: the systemic velocity 
corrected to the heliocentric reference frame V s =3587 + 9 km s _1 , 
¥0 =128° + 6°, M =2.7+0.3 x 10 9 M Q , i = 38 + 5°, and 
A =238 + 15 pc. The derived kinematical centre is coincident with 
the adopted position for the nucleus (peak of the continuum emis- 
sion) within the uncertainties, with X = -17+7 pc and 7 =3+3 pc, 
considering that our spatial resolution is (FWHM) ^0715 (35 pc). 
The uncertainty in the parameters were obtained directly from the 
application of the Levenberg-Marquardt method, taking into ac- 
count the uncertainty in the velocity. 

The parameters derived from the stellar kinematics can be 
compared with those from previous works. The systemic velocity 



10 Riff el and Storchi-Bergmann 




-17.0 



-17.2 



-17.5 



-17.7 



-18.0 



-18.2 



-18.4 



-18.7 



-18.9 



Figure 10. Velocity slices along the H2 i2.1218//m emission-line profile. Each panel shows the intensity distribution at the velocity shown in the top-left 
corner of the panel, integrated within a velocity bin of ~60 kms -1 . The cross marks the position of the nucleus, the white line shows the galaxy major axis, 
the arrows identify regions where we identify an asymmetry in the flux distribution relative to the major axis, and the green contours overlaid to some panels 
are from the 3.6 cm radio image. 



derived here is abo ut 40 kms -1 smaller than the one obtained by 
iBower et"aD (1 1995b from optical emi s sion li nes and in agreement 
with the one derived by iKnop et all d200lh using near-IR emis- 
sion lines. We consider our measurement of V s more robust than 
those from previous works as our observations have higher spatial 
resolution and V s was obtained from the stellar kinematics, which 
is a better tracer of the gravitational potential of the galaxy than 
the gas. The orientation of the line of nodes ¥0 =128° i s in re a- 
sonable agreement with the one obtained bv lKnop et all (1200 ll) - 
¥0 =120°. The scale length A and the bulge mass M are similar to 
those obtained for the ce ntral region of other Seyfert galaxies us- 
ing similar modeling (e.g. lBarbosa et al.l20 06). It should be noticed 
that M is tightly coupled with the inclination (/) of the disc, since 
V 2 oc Msin(/). Thus the uncertainties in M and i can be even higher 
than those quoted above due to this degeneracy. Nevertheless, the 
value obtained for the inclinati on of the disc is in r easonable agree- 
ment with the one obtained bv lBower et al.l i 19951) . i = 45°. 

In figure QT] we show the best fit model for the stellar velocity 
field (top panel) and the residual map, obtained from the subtraction 
of the model from the observed velocities. A comparison between 
the modeled and the observed velocities is also shown in the top-left 
panel of Fig. [12] where we present one-dimensional cuts for the ve- 
locities along the line of nodes, averaged within a pseudo-slit (J/25 



wide. The bottom panel of figure [TT] shows that the residuals be- 
tween the observed and modeled stellar velocity fields are usually 
smaller than 20 km s _1 (^ 15% of the maximum velocity amplitude 
of the rotation curve), but present a systematic behaviour. These 
systematic residuals - blueshifts to the south-west of the major axis 
and redshifts to the north-east in the inner arc second radius - seem 
to correlate with the oval seen in the HST F606W optical image, 
as shown by the corresponding superimposed contours. A signa- 
ture of the effect of this oval is already present in Fig. [2] where the 
"zero-velocity curve" shown shows an S shape, a known signature 
of the presence of a non-axisymetric structure in the gravitational 
potential, which is not present in the model stellar velocity field of 
Fig.fTTlwhich shows straight isoveloci ty lines for velocities cl ose to 
zero (along the kinematic minor axis uEmsellem et all d2006h have 
found a similar S- shaped structure in the stellar velocity field of the 
Seyfert galaxy NGC 1068, attributed to the presence of a nuclear 
bar. In the case of Mrk 1066, we attribute the distortions to the nu- 
clear oval seen in the the HST image from which spiral arms seem 
to originate along PA=128° (see Paper I). 

This oval is also seen in the stellar velocity dispersion map 
(top left panel of Fig. [2]), in the form of a partial circumnuclear ring 
of low values with cr* & 50 km s _1 , at a distance of « 1" from the 
nucleus, immersed in a background with values of cr* « 90 km s _1 . 
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Figure 11. Rotating disc model for the stellar velocity field of Mrk 1066 
(top panel) and residual map between the observed and modeled velocities 
(bottom panel). The dashed line shows the orientation of the line of nodes 
OFo =128°), the cross marks the position of the nucleus and the contours 
are from the HST optical image. 



This partial ring is also observed in the one-dimensional cut of the 
<t* map shown in the top-right panel of Fig. [12] Similar circumnu- 
clear rings of low velocity dispersions have been reported for sev- 
eral o ther Seyfert galaxies fearbosa et al]l2006l : liiffel et alJl2008L 

2009) , interpreted as colder regions with more recent star forma- 
tion than the underlying bulge. This interpretation is su pported by 
a recent study of the stellar population of this region (Riffe l etal.1 

2010) , which shows a dominance of intermediate age (10 8 - 10 9 yr) 
stars. 

The /* 4 * map (bottom-right panel of Fig. shows a ring of 
high values co- spatial with the ring of low <x* values. Such /z 4 * en- 
hancement (due to a velocity distribution which is narrower than a 
Gaussian for low velocities) can be attributed to the contribution of 
the youngest stars (the stars actually span a range of intermediate 
ages) in the ring, which have a colder kinematics than the surround- 
ing bulge stars. 

The mass of the super-massive black hole in the centre of 
Mrk 1066 can be estimated from the bulge stellar velocity dis- 
persion using \og(M BH /M Q ) = a + /?log(o%/cr q), where a = 
8.13 + 0.06,/? = 4.02 + 0.32 and cr = 200 km s" 1 dTremaine et aD 
120021) . We adopt cr* « 90 kms -1 as representative of the bulge, 
a value in good agreem ent with previous optical measurements 
dNelson & W hittle 1995), and obtain M BH = 5.4^ x 10 6 M Q . 



4.2 Gas kinematics 

Although the gaseous velocity fields of Fig.|4]present also a similar 
rotation pattern to that observed for the stars, they are more dis- 
turbed, indicating that additional kinematic components - not dom- 
inated by the gravitational potential of the bulge - are present. Dif- 
ferences between the gas and stellar kinematics are also observed in 
the one-dimensional cuts along the major axis of the galaxy, shown 
in the top panels of Fig.[l2] To south-east of the nucleus the gaseous 
velocities are redshifted by ^30 kms -1 relative to the stellar ones, 
while at «075 north-west they are redshifted by ^40 km s _1 . Some 
differences between the [Fen], H 2 and Pa/? kinematics are also ob- 
served, in particular between 1" and 1"5 north-west from the nu- 
cleus, where the H 2 and Pa/? velocities are similar to the stellar and 
the [Fen] velocities are blueshifted by ^30 km s _1 . A distinct kine- 
matics for the H 2 , [Fen] and Pa/? is also supported by the cr maps 
(Fig. [5]) and one-dimensional cr cuts (top-right panel of Fig. Q2}, 
which show that, at most locations, the H 2 presents the smallest cr 
values, followed by Pa/?, with [Fe 11] presenting the highest cr val- 
ues. Differences between the kinematics of the ionized and molec- 
ular gas are also observed in the h 3 and h$ maps (Figs. [6] and |7] 
respectively), as well as in the channel maps (Figs. [8] [9] 013, which 
shows [Fen] emission at highest blue- and redshifts than H 2 , indi- 
cating that the H 2 emission traces less disturbed gas than the [Fe 11] 
emission. 

The above c onclusion is supported by other recent works. 
Hick s et all d2009h have studied the H 2 kinematics from the inner 
^100 pc of a sample of 9 Seyfert galaxies using SINFONI IFU 
at the ESO Very Large Telescope (VLT) and concluded that it is 
dominated by rotation in a disc with typical radius of ^30 pc and a 
comparable height. In previous studies by our group we have also 
found that while the [Fen] emitting gas has important kinematic 
components attributed to gas extending to high galactic latitudes 
and in interaction with a radio jet, the H 2 kinematics is dominated 
by rotation in the galaxy disc jRodnguez-Ardila. Riffel & Pastorizal 
| 2005 |: lRiffel et alJ200d E008. 2009: IStorchi-Bergmann et alll999L 
l2009h . 

We can compare our results w ith previously publ ished near- 
IR dKnop et all 1200 lh and optical feower et al J 1 19951) gas kine- 
matics for Mrk 1066 using long-slit spectroscopy. iBower et"aD 
dl995h present measurements of the centroid velocity and FWHM 
of the Her, H/?, [O m] ^5007 A, [Oi] ,16300 A, [N 11] ^6584 A, 
[S 11] /16724 A emission lines along PA=134°, obtained under a see- 
ing of «079. All emission lines, with exception of the [O m], present 
similar velocities which are consistent with our near-IR measure- 
ments. Nevertheless, our velocity field reveals much more details, 
which we attribute to our better spatial resolution. The [O m] emis- 
sion is blueshifted relative to the other emission lines to north- 
west of the nucleus and redshifted to south-east. The presence of 
this kinematic component in our data is suggested by our h 3 maps 
for [Fe 11] and Pa/? (bottom-left and top-right panels of Fig. [6} and 
one-dimensional cuts (bottom-left panel of Fig. [12]), which show 
positive values (red wings) to south-east of the nucleus and nega- 
tive values to north-w est (blue wings). Th e presence of wings have 
also been observed bv lKnop et al. 1 12001) for the near-IR em ission 
lines along PA=135° at a seeing of a/7. iKnop et all (1200 ll) sug- 
gested in addition that t hese wings a r e asso ciated with the [O m] 
kinematic component of IBower et al.l (1 19951) . In particular, at «1" 
north-west of the nucleus the [Fen] emission line is blueshifted 
by ^30 kms -1 relative to the other near-IR emission lines, sug- 
gesting that the [Fe 11] emitting gas has the kinematical component 
observed in [O in]. In addition, both present higher velocity disper- 
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Figure 12. One-dimensional cuts for the stellar and gaseous kinematics obtained along PA= 128° (orientation of the line of nodes) within a pseudo-slit 
0725 wide. Top: Radial velocities (left) and velocity dispersions (right). Bottom: (left) and h\ (right) Gauss -Hermite moments. The model is shown as a 
continuous line in the top-left panel. Typical uncertainties for <x are <10 kms -1 and for and h\ are « 0.03. 



sions when compared to other optical and near-IR emission lines 
with peak values observed at «1" north-west of the nucleus, close 
to a radio hot spot. These results support the conclusion that at least 
part of the [Fe n] emitting gas has the same kinematics of the [O m] 
emitting gas. 

As our goal here is to look for non-circular motions in the gas, 
we have subtracted the stellar velocity field (which is dominated 
by circular rotation) from the gaseous ones, considering also that 
we could not fit a circular model to the gas kinematics, which is 
too much disturbed. Figure [13] shows the resulting maps for the the 
[Pn] (top-left panel), [Fen] (top-right), Pa/? (bottom-left) and H 2 
(bottom-right). Although these residual maps should be considered 
with caution, as the gas respon se to the gravitation al potential dif- 
fers from that of the stars (e.g. lCombes et al.ll 19958) , we have used 
them only as guides, together with the gas centroid velocity fields 
and channel maps (Figs. [4] [12] [8] [9] and [TO]) in order to look for 
non-circular velocity components. Considering all these maps, we 
tentatively identify the presence of two structures in the gas (be- 
sides the gas in similar rotation to that of the stars): (i) a compact 
disc with excess velocity relative to that of the stars of ^40 km s _1 
within ^074 from the nucleus, more clearly observed in the H 2 
residual map (bottom right panel of Fig. [13]); (ii) an outflow identi- 
fied by excess blueshifts of ^50 kms -1 at 173 to north-west of the 



nucleus, close to the edge of the radio jet, more clearly observed 
in the [Fe n] and [P n] residual maps (top left and right panels of 
Fig. [13]). A similar structure in redshift is observed at «1" south- 
east of the nucleus. We now discuss the observational signatures of 
these two components - the compact disc and the outflow. 



4. 2. 1 The compact disc 

This structure is observed not only in the residual map, but also 
directly in the H 2 velocity field (bottom-right panel of Fig. [4]), 
which shows an abrupt increase in velocity between the nucleus 
and « 073 (« 70 pc), while the stellar rotation shows a much shal- 
lower increase in the rotation velocity, what can also be seen in the 
one-dimensional cuts of Figure [12] The presence of this disc in H 2 
is further supported by the low velocity dispersions (see Figs.[5]and 
[T2lh with values smaller than the stellar ones, observed closer to the 
nucleus in H 2 than in the other emission lines. 

Recent IFU observations of other Seyfert galaxies by other 
authors have also revealed the presence of circum nuclear H 2 discs 
drlicks et alJl2009h . Following Ifficks etal.1 d2009h . we estimate the 
dynamical mass within the compact disc radius (70 pc) from its 
kinematics as: 
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Figure 13. Residual maps for the gas kinematics - differences between the gaseous and stellar velocity fields for [P 11] (top-left panel), [Fe 11] (top-right panel), 
Pa/5 (bottom-left panel) and H2 (bottom-right panel). The thick contours overlaid to the [Fen] map are from the radio image, while those overlaid to the H2 
map are from the HST optical continuum. The continuous line marks the line of nodes. 



Mdyn M VklpL, (5) 

where V rot = V bs/sin(/) is the rotational velocity, V bs is the ob- 
served velocity, G is the gravitational constant. From Fig. [12] we 
observe that at 073 from the nucleus, the H 2 V bs ~65 kms -1 (av- 
erage of the values observed to north-west and to south-east of the 
nucleus) and cr « 60km s" 1 , resulting in M dyn « 3.6 X 10 8 M . As - 
suming a typical gas mass fraction of 10% fsee lfficks etall (120091) 
for a discussion about the observed mass fraction), we estimate that 
the gas mass in the disc is M gas « 3.6 x 10 7 M . This m ass is about 
one or der of magnitude larger than those obtained by iHicks et all 
(2009) for a sample of 9 active galaxies, what can be due to the fact 
that the radius of the nuclear disc in Mr k 1066 is at least two times 
larger than those in IHicks et al.1 d2QQ9h . The mass derived above 
probably consists mostly of cold molecular gas since, as pointed 
out in Paper I, the mass of hot molecular gas (which emits in the 
near-IR) is about 10 4 times smaller than the dynamical mass ob- 



tained here. Radio observations indeed show the presence of large 
amounts of cold molecular gas with masses o f 10 7 -10 9 M in the 
inner hundred parsecs o f active galaxies ( e.g. iGarcfa-Burillo et aD 
l2005l : lBoone et al.ll200llkrips et al.ll2007h . 

Feeding of the compact disc 

As pointed out in Section [331 the H 2 channel maps (Fig. [TUt 
are not symmetric with respect to the major axis in the blueshifted 
channels centred at -104 kms -1 and -42 kms -1 , where the flux 
distribution curves towards the north, and in the redshifted channels 
centred at 81 km s -1 and 142 km s -1 where there is a similar curva- 
ture to the south (see the arrows in Figure [10]). In order to interpret 
these features, we now consider that the major axis of the galaxy is 
at PA=128°, running from the bottom left to the upper right in the 
panels of Fig.[l0] as shown by the white line. Under the assumption 
that the spiral arms (outlined by the contours of the HST image su- 
perposed in Fig. [2) are trailing, we conclude that the north-east is 
the far and the south-west is the near side of the galaxy. If the H 2 
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Structure Map 




Figure 14. Structure map of the inner region of Mrk 1066 obtained from 
the HST image. The arrows have been drawn at the locations they appear in 
Fig. [TO] and approximately trace nuclear spiral arms. 



emitting gas is in the galaxy plane, as suggested by its low veloc- 
ity dispersion (bottom right panel of Fig.[5] and upper right panel 
in FigIT2l we suggest that the curves in the light distributions in 
the channels maps centred at -104kms _1 and -42kms _1 to the 
north-east, as well as in those centred at 81 kins" 1 and 142kms _1 
to the south-west (see arrows) are due to inflows along spiral arms. 
We note that the same structures - resembling spiral arms - are 
observed also in the H 2 velocity dispersion map (Fig. [5]). 

Morphological nuclear spiral arms can also be seen 
in the optical HST F606W image of Mrk 1066 (see 
RiffeL Storchi-Bergmann & Nagar] l2010h . A structure map, 
obtain ed from this image using the technique of Pogge & Martinil 
(2002) is shown in Fig.PHl The arrows in this map approximately 
trace dark spiral arms and correspond to the same location we have 
drawn them in the H2 channel maps. In the structure map, the dark 
spiral arms are due to dust. It is expected that the H 2 is associated 
to dusty regions, where it is protected from strong radiation which 
may destroy the H 2 molecule. This has been shown to be the case 
for NGC 1097 bv lDavies etaD (120091) . who found velocity residu- 
als in the H 2 kinematics which were ass ociated with dusty spirals 
(Prieto, Maciejewski & Reunanen l2005h . interpreted as inflows 
towards the nucleu s. Inflows were als o observed along these spiral s 
in ionized gas by iFathi et al ] d2006h . Is imdes Lopes et al.1 b007h . 
using HST optical images for a matched sample of galaxies with 
and without AGNs, found that, for the early-type subsample, all 
AGN hosts have circumnuclear structures of gas and dust, while 
this is observed in only 26% of the i nactive galaxies. In many 
cases, the dusty structures observed bv lSimoes Lopes et al.1 d2007h 
resemble spiral arms, suggesting that they are the feeding channels 
of the AGN, as confirmed for a few cases in which the kinematics 
has been obtained. 

Inflows alon g nuclear spirals have been predicted in numeri- 
cal simulations by Maci eiewskil j2004al lbL who demonstrated that, 
if a central SMBH is present, spiral shocks can extend all the way 
to the vicinity of the SMBH and generate gas inflows consistent 



with t he accretion rates inferred in local AGN. iMartini & Poggel 
(1999) in their study of nuclear spirals in HST images have in- 
deed concluded that they are not self-gravitating, being consistent 
with shocks in nuclear gas discs. Kinematical spiral arms have 
also been observed in the central region of other Seyfert galax- 
ies (IFathi et all 120061: IStorchi-Bergmann et al 1 120071: iRifFel et al ' 
2008 : ISanchez et al.ll2009l : Ivan de Yen & Fathill2009l : iDavies et al 



2( 
2( 



2009). The modeling of the gaseous velocity field by a rotating disc 
plus spiral shocks shows that the gas is slowly streaming towards 
the nucleus with a typical mass inflow rate of lO^-lO 1 M /y r 



(IDavies et al.l2009l : lvan de Yen & Fathil2009l : ISanchez et al.l200 ( 
Thus, the nuclear spiral arms seen in the H 2 kinematics in Mrk 1066 
could also be the feeding channels of its SMBH, or at least of the 
compact disc discussed above. 

4.2.2 The outflowing gas 

The top-right panel of Fig. [T3l shows blueshifts in the [Fen] resid- 
ual map of up to -80km s _1 at ~ 1" north-west of the nucleus, at 
the top edge of the radio jet. Similar residuals are observed for the 
[P n], and somewhat smaller for Pa/? and H 2 (Fig. [13]). To the south- 
east the residual maps present similar redshifts in a region near the 
bottom edge of the radio jet. We interpret these residuals as due to 
an outflow and in particular from regions of interaction between the 
radio jet and the ambient gas. This interpretation is supported by (i) 
the residuals being observed close to the ends of the radio jet, as 
mentioned above; (ii) the fact that the [Fen] emitting gas presents 
increased velocity dispersion values in regions surrounding the ra- 
dio structure. Such cr enha ncements are expected when gas i s dis- 
turbed by a radio jet (e.g. iDopita & Sutherland! Il995l [l999); (iii) 
The h 3 map (Fig. [6]) for the [Fe 11] presents its most negative values, 
of « -0.3, in a region co- spatial with the blueshifts observed in the 
residual map; (iv) The velocity channel maps of [Fen] are observed 
up to higher blueshifts and redshfits than those for Pa/? and H 2 . The 
flux distributions at these high velocities are well correlated with 
the radio structure. This behaviour is similar to that observed for 
the [O m] emission (iBower et~aDll995h . interpreted as originated in 
outflowing gas in a bi-conical structure with axis approximately co- 
incident with the orientation of radio jet, while the lower excitation 
optical line emission is attributed to gas located in the plane of the 
galaxy. 

The outflowing component is not more obvious in the near-IR 
because of the orientation of the outflow, which is close to the major 
axis of the galaxy, where the gas kinematics is dominated by rota- 
tion in the plane (as observed in the stellar velocity field of Fig. 13 
and in the compact disc (Sec l4.2.iT l. The blueshifts described above 
represent a kinematic component which is superimposed on the 
galaxy rotation and on the compact disc component. The blueshifts 
at ~ 1" north-west and redshifts at similar distances to the south- 
east can be interpreted as originating in a bi-conical outflow which 
shares the [O m] kinematics but is outshined by emission from the 
galaxy disc and from the compact nuclear disc in the inner regions. 

The bi-conical outflow is better observed in the [Fe 11] chan- 
nel maps (Fig. [8]), which show the highest blueshifts (from -300 to 
-500 km s _1 ) to north-west of the nucleus and the highest redshifts 
observed to both sides of the nucleus. Comparing the channel maps 
with the [Fe 11] velocity field (Fig. [4]) we conclude that the north- 
west blueshifts and south-east redshifts are due to outflows, as the 
galaxy rotation is in the opposite direction, showing blueshifts to 
the south-east and redshifts to the north-west. At intermediate ve- 
locities (80-220 kms -1 ) both redshifts and blueshifts are observed 
to both sides of the nucleus confirming the presence of two velocity 
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components: rotation and outflows. A similar behaviour is observed 
for the Pa/? channel maps. Lower velocity channels are consistent 
with emission from gas in the plane of the galaxy. The H2 chan- 
nel maps are dominated by these lower velocities, supporting the 
interpretation that the H 2 is located in the galaxy plane. 

Under the assumption that part of the [Fen] emitting gas 
shares the [O m] kinematics, we can estimate the opening angle of 
the cone directly from Fig.[l3]as being a « 20°. As the radio jet is 
oriented almost along the galaxy major axis, we can constrain the 
angle of the bi-cone axis relative to the plane of the sky as being 
close but somewhat larger than zero (6 > 0°), since the north-west 
side of the bi-cone must be oriented towards us (although appar- 
ently by a small angle) in order to produce excess blueshifts there 
while the south-east side of the bi-cone is directed away from us in 
order to produce the observed redshifts. 

We now estimate the mass outflow rate across a circular cross 
section along the bicone located at 1"3 from the nucleus, which 
has a radius r ^0723 « 54 pc for a bi-cone opening angle of 20°, 
corresponding to an area of A « 8.7 x 10 40 cm 2 . The component 
of the velocity of the outflowing gas along the bi-cone axis (v out ) is 
related to the observed velocity (v obs ) by v out = v obs /sin#, where 6 
is the angle between the bi-cone axis and the plane of the sky, and 
thus the mass outflow rate can be obtained from 



M out = 



2m p N e v ohs fA 
sin 6 



(6) 



where the factor 2 is included in order to account for the outflows 
to both sides of the nucleus. Assuming N e = 500 cm -3 , / = 0.01 
and v out = 50 kms -1 , directly form Fig. |T3j we obtain M out « lx 
10- 2 /sin6>M o yr" 1 . 

In order to estimate the value of 6, we consider the observation 
that the [O m] image presents a collimated structure t o north-west 
of the nucleus and a fainter emission to the south-east teower et alj 
1995), suggesting that the bi-conical outflow makes a small angle 
with the plane of the sky, with the south-east side apparently ob- 
scured by the galaxy plane. A large angle would produce a less 
collimated emission. 

In order to quantify the mass outflow rate, we assume 6 = 
10°. In this case the adopted geometry implies that the far wall 
of the north-west cone and the near wall of the south-east cone 
are exactly in the plane of the sky. With this value of 6 we obtain 
M out * 6xl0- 2 Moyr" 1 . 

The mass outflow rate obtained for Mrk 1066 can be 
compared with those obtained for other active galaxies. 
Veill eux. Cecil & Bland-Hawthorn! (2005) report values of M out « 
0.1 - 10 M yr _1 for warm ionized gas ou tflow from luminous 
active galaxies, while iBarbosa et al.l (120091) report M out & (1 - 
50) X 10 -3 M© yr -1 from optical IFU observation of a sample of six 
Seyfert galaxies assuming and electron density of N e = 100 cm -3 . 
IStorchi-Bergmann et al.1 f2010) found M out « 2 M yr _1 from the 
nucleus of NGC4151 us ing similar NIFS observation s of its cir- 
cumnuclear region, while Crenshaw & Kraemerl (120070 obtained an 
outflow rate « 6 times smaller using blueshi fted absorp tion lines 
in the UV spectrum of NGC4151. In iRiffel et aD d2009h we es- 
timated M out « 5 x 10" 2 M yr _1 for the Seyfert 2 nucleus of 
NGC 7582 for the ionized gas using K-band IFU spectroscopy. The 
mass outflow rate estimated for Mrk 1066 is thus within the range 
of observed values for active g alaxies, being 2 times smaller than 
the lo west ones obtained by IVeilleux. Cecil & Bland-Hawthornl 
2005 ), « 30 times smaller than the one obtained for NG C 4151 



Storchi-Bergmann et al . 2010; Crenshaw & Kraemer 2007), about 



10 times larger than the values obtained bv lBarbosa et al. I d2009h 
and similar to the outflow rate of NGC 7582 (Riffel et al. 2(X)i). 



4.2.3 Feeding vs feedback in Mrk 1 066 

iKnop et all (1200 lh - using long- slit spectroscopy - proposed that 
most of the near-IR line emission in Mrk 1066, as well as the optical 
lines of [N n], [S n], [0 1] and H 1, originate in a rotating disc in the 
plane of the galaxy with a major axis oriented along PA^ 120°, 
while the wings of the near-IR emission lines and the optical [O m] 
emission originate from the bi-conical outflow, as illustrated in their 
Fig. 8. 

Our observations have 2D cove rage, as well as higher spectral 
and spatial resolutions than those of IKnop et al.l d200lh and can be 
used to better constrain the orientation of the cone and major axis 
of th e disc and we can thus further constrain the scenario proposed 
bv IKnop et al.1 J200lh. In Fig. Q3] we present a cartoon illustrat- 
ing our proposed physical scenario for the circ umnuclear region o f 
Mrk 1066. This cartoon is inspired on Fig. 8 o f Knopet al. (2001), 
in which the high ionization emission lines - such as the optical 
[0 111] lines - originate in outflowing gas along the bi-cone oriented 
along the radio axis (PA^ 135/315°) and with an opening angle of 
^20°. 

In agreement with lKnope t al. (2001), we support that the low- 
excitation line emission (from H 2 , Hi, [Fen], [P11], [Nil], [Oi], 
[S 11]) is dominated by gas rotating in the galaxy disc, with similar 
kinematics to that of the stars: the line of nodes is oriented along 
PA « 128° and the disc inclination is i « 38° relative to the plane 
of the sky. But we can also identify in the near-IR emission lines a 
kinematic component originating in the bi-conical outflow (which 
dominates the [O m] emission), as evidenced by the blueshifts to 
north-west and redshifts to south-east of the nucleus in Fig. \\3\ 
The contribution of this component to the gas kinematics increases 
in importance from H 2 — >Pa/?— >[Fen]. This can also be observed in 
the channel maps: while the presence of outflows combined with 
rotation is clear in the [Fe 11] and Pa/? channel maps, the H 2 chan- 
nel maps seem to show only rotation. These kinematics suggest 
that part of the ionized gas extends to high galactic latitudes, where 
it contributes to the bi-conical outflow, while the molecular gas is 
mostly restricted to the galaxy plane. 

This picture - in which the H 2 originates mostly from the 
galaxy plane, while the ionized gas is observed also in outflows 
extending to high galactic latitudes - is in good agreeme nt with re- 
sults we have obtained for other Seyfert galaxies (e.g. IRiffel et all 
2006, 2008, 2009; Storchi -Bergmann et al . 1999, 2009, 2010). 

Another interesting new result is the observation of a compact 
nuclear disc, with radius ^70 pc, which is rotating faster than the 
underlying stellar disc. We interpret this nuclear disc as a gaseous 
structure which will soon give origin to the formation of new stars. 
We are probably witnessing the formation of a nuclear stellar disc, 
such as those observed as low stella r velocity dispersio n structures 
surrounding active galactic nuclei ( Davies et al.l 120071) or as high 
surface brightness structures in previous HST imaging studies (e.g. 
Simoe s Lopes et al ] |2007l : IPizzella~e t al. 2002). This compact disc 
is apparently being fed by molecular gas flowing along nuclear spi- 
rals, as evidenced by the blueshifts observed at the far side of the 
galaxy and redshifts in the near side (Fig.[l0]> si milarly to what w e 
have observed in the Seyfert galaxy NGC 4051 (IRiffel et al.ll2008h . 

Finally, we now compare the outflow mass rate, obtained in 
Sec l4.2.21 with the accretion rate necessary to power the AGN at 
the nucleus of Mrk 1066. 
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Figure 15. Physical scenario for the circumnuclear region of Mrk 1066 (iKnop et al1l200lh . Most of the near-IR line emission originates in a rotating disc with 
major axis coincident with that of the stellar disc (white dashed line along PA=128°). The [Om] line emission originates mostly from gas in the bi-conical 
outflow with redshifts to south-east and blueshifts to north-west of the nucleus, opening angle « 20° and oriented along the PA of the radio jet (arrows along 
PA^ 135/3 15°). The outflowing component observed in the near-IR emission lines also originates in the bi-conical outflow. The arrows shown along the minor 
axis represent the nuclear spiral arms observed in Fig. [13] 



The mass accretion rate to feed the active nucleus can be ob- 
tained from 



where Lb i is the nuclear bolometric luminosity, 77 is the efficiency 
of conversion of the rest mass energy of the accreted material 
into radiation and c is the light speed. The bolometric luminos- 
ity can be approximated by L bo i » 100L na , where Ln a is the Ha 
nuclear luminosity (e.g. iHol [l999l : iHo et all l200lh . The Bry nu- 
clear flux measured for an circular aperture with 0725 radius is 
F Brr ~ 4.2 x 10~ 15 erg s _1 cm -2 . In Paper I, we obtained a red- 
dening of E(B - V) « 1 for the nucleus of Mrk 1066 using 
the Bry/Pa/3 line ra tio. Correcting the Bry flux for t his redden- 
ing using the law of Cardelli, Clayton & Mathis (1989), we obtain 
F Brr ~ 5.8 x 10~ 15 erg s _1 cm -2 . Assuming a temperature T = 10 4 K 
and an electron density n e = 10 2 cm -3 , the ratio between Ha 
and B ry is predicted to be Fu a /F Bri/ = 103 dOsterbrock & Ferlandl 
l2006h . Thus L Ua * 1.7 X 10 41 erg s" 1 for a distance to Mrk 1066 of 
d = 48.6 Mpc and the nuclear bolometric luminosity is estimated 
to be L bo i ~ 1.7 X 10 43 erg s _1 . Assuming 77 ^ 0.1, which is a typical 
value for a "standard" geometrically thin, optically thick accretion 
disc (e.g. lFrank. King &~R aine 2002), we obtain a mass inflow rate 
of m * 3 x 10" 3 Moyr 1 . 

The mass outflow rate (M out ) in the NLR is about one order 
of magnitude larger than m, a ratio comparable to those observed 



for other Seyfert galaxie s (Barbosa et al. 2009; Riffel et al. 2009; 
IStorchi-Bergmann et al.ll201ok which indicates that most of the 
outflowing gas in the NLR of active galaxies does not originate 
in the AGN but in the surrounding interstellar medium from the 
galaxy plane, which is pushed away by the nuclear outflow. 



5 SUMMARY AND CONCLUSIONS 

We used integral field J and Ki bands spectroscopy from the inner 
^350 pc radius of the Seyfert galaxy Mrk 1066, obtained with the 
Gemini NIFS at a spatial resolution of ~35 pc and spectral resolu- 
tion of ~40kms _1 , to map the gaseous and stellar kinematics. Our 
main conclusions are: 

• The stellar velocity field is dominated by circular rotation in 
the plane of the galaxy, but shows an S- shape distortion along the 
galaxy minor axis which we attribute to the presence of an oval 
structure in the galactic potential. 

• The oval structure is also delineated by a partial ring of low 
velocity dispersion (cr* ^50kms _1 ) at ^230pc from the nucleus, 
which is due to intermediate- age stars immers ed in the old bulg e, 
as concluded from stellar population synthesis teiffel et aDl2010h . 

• The bulge velocity dispersion is cr* « 90kms _1 , implying in 
a black hole mass of M BH ~ 5.4 x 10 6 M Q . 

• The gaseous kinematics shows circular rotation similar to that 
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of the stars with two additional kinematic components: (i) a com- 
pact nuclear disc, more clearly observed in H 2 emission, with radius 
r = 70 pc and gas mass of ~ 10 7 M Q ; (ii) outflows associated with 
the radio jet, presenting a mass outflow rate of « 6 x 10 _2 M o yr _1 , 
more clearly observed in [Fe n] emission. 

• The compact H 2 nuclear disc presents the smallest gas cr val- 
ues (< 70 km s _1 ), consistent with an origin in the galaxy plane and 
is apparently being fed via nuclear spiral arms, observed in molec- 
ular gas at similar velocity dispersion to that of the compact disc. 
This disc may not only be the feeding source of the AGN but prob- 
ably will be consumed in star formation, leading to the formation 
of a nuclear stellar disc, a structure frequently observed in galaxy 
bulges. 

• The [Fen] emission presents the largest cr values (up to 
150 km s _1 ) and the highest blue- and redshifts, of up to 500 km s _1 , 
while the highest stellar rotation velocity is only ^130 km s _1 . The 
high velocities and velocity dispersions are found at the tip and 
around the radio jet, which is oriented approximately along the 
galaxy major axis, with the north-west side slightly tilted towards 
us. The outflow rate is about one order of magnitude higher than 
the accretion rate, indicating that most of the outflow consists of 
entrained gas from the galaxy plane by a nuclear outflow. 
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